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O , Abstract 

in ■ 

Q\ ■ We study the scattering of hard external particles in a heat bath 

in a real-time formalism for finite temperature QED. We investigate 
Q_i| the distribution of the 4-momentum difference of initial and final hard 

O ■ particles in a fully covariant manner when the scale of the process, Q, 

is much larger than the temperature, T. Our computations are valid 
for all T subject to this constraint. We exponentiate the leading infra- 
red term at one-loop order through a resummation of soft (thermal) 
^ ■ photon emissions and absorptions. For T > 0, we find that tensor 

structures arise which are not present at T = 0. These carry thermal 
signatures. As a result, external particles can serve as thermometers 
introduced into the heat bath. We investigate the phase space origin 
of log(Q/m) and \og(Q/T) terms. 
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1 Introduction 



Recently there has been a lot of phenomenological interest in external par- 
ticles interacting inside a heat-bath. This is due to the realisation that hard 
jets formed in the early stages of a heavy-ion collision can interact with a 
quark-gluon plasma formed later, and perhaps carry away information on 
this plasma Jj], 0. However, such processes are often infra-red sensitive. 
Phenomenological studies have used arbitrary infra-red cutoffs to obtain fi- 
nite results. It is legitimate to ask whether one can remove this arbitrariness 
and obtain stable finite results for such infra-red sensitive processes. 

In this paper we consider the simpler, but related, problem of hard scat- 
terings of some external particles introduced within a heat-bath at a temper- 
ature T. We choose two of the simplest processes in QED — the scattering 
of an electron with an off-shell photon and the annihilation of an electron- 
positron pair into an off-shell photon. Due to quantum and thermal fluctu- 
ations, extra soft particles can also be involved in each of these scatterings. 

If p m and p ont are the summed momenta of all the hard initial and final 
state particles respectively, then one can construct the difference 

p»=p»-pT- (i-i) 

We study the distribution of P through the measurement of a suitable semi- 
inclusive cross section. It is non-trivial because of the extra soft particles. 
When P 2 <C Q 2 (Q 2 is the scale of the hard process), the Lorentz-invariant 
cross section factorises as 

d 4 a d 4 V 



dp4 " (1-2) 

where o"o is the inclusive cross section for the process. Throughout this paper 
we assume such a factorisation, and study the distribution of P M (the second 
factor on the right hand side) through its cumulants, i.e., the connected parts 
of the moments. 

Naively, it might seem that the heat-bath spoils the Lorentz invariance 
of various cross sections. This is incorrect. The bath only selects a frame; 
Lorentz invariance cannot be lost due to the existence of a special frame. 
Cross sections which are Lorentz invariant at T = remain so for T > 0. We 
carry out a fully covariant treatment of the resummation and show that the 
difference between T = and T > physics lies in some new asymmetries and 
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angular correlations. These are proportional to powers of the temperature. 
As a result, scatterings of external particles can be used as thermometers for 
the heat-bath. This is the main result of this paper. 

It is useful to recall the status of infra-red sensitive cross sections at T = 0. 
Exclusive scattering cross sections can be infra-red singular. The analysis due 
to Bloch and Nordsieck |3j showed that such measurements are impossible to 
make with any real detector. They argued that a real detector, with a finite 
energy resolution, measures, at best, some semi- inclusive cross section. They 
computed one such cross section and showed that it is indeed non-singular. 
In the limit where the detector becomes ideal, i.e., the resolution becomes 
infinitely sharp, this cross section remains non-singular. 

A stronger theorem is also available for T = field theories; it goes by the 
name of the Kinoshita-Lee-Nauenberg theorem The theorem states that 
inclusive cross sections are infra-red finite, order by order in perturbation the- 
ory. Infra-red singularities in real and virtual diagrams cancel in each order. 
Semi-inclusive cross sections are also infra-red finite. However, incomplete 
cancellations between real and virtual parts show up as large logarithms. It 
is possible to sum these to all orders in perturbation theory. 

At finite temperatures the infra-red singularities could be stronger, since 
every squared matrix element comes with a thermal distribution function 
attached to each external leg. Each external boson then adds one power 
to the infra-red divergence. A T > version of the KLN theorem has not 
been proved. Computations have been performed in several models to check 
whether or not infra-red divergences cancel in suitably defined inclusive or 
semi-inclusive processes || |J. 

An important recent result is due to Landshoff and Taylor 0. In a 
leading order perturbative computation of the decay width for a scalar going 
to two scalars, they found exact cancellation of the infra-red singularity. 
Weldon introduced a crucial simplification in finite temperature QED ||. By 
eikonalising the soft part of the one-loop order cross section, he was able to 
resum soft photon effects to all orders in perturbation theory and obtain the 
distribution dV/dPo. This is an extension of the Bloch-Nordsieck theorem 
to finite temperature. 

This proof is a crucial input to our work. We extend it first to the case 
of interest to us, i.e., the full momentum resummation. We construct a 
Lorentz invariant characteristic function for the distribution of P, by expo- 
nentiating the first order perturbative correction to the hard process. Then 
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we extract the cumulants of this distribution and find terms of the form 
log(Q/m). These large logarithms come from parts of phase space where the 
extra photons are nearly collinear with the fermion momenta, exactly as at 
T = 0. These large terms allow us to identify phenomenologically interesting 
results from such soft photon resummations at finite temperature. We also 
find terms in log(Q/T). These arise from terms left over after cancellation 
of the infra-red divergence. 

The plan of this paper is the following. In Section |2] we shall perform 
the resummation and show that Weldon's proof of infra-red finiteness |§ can 
be generalised in a Lorentz invariant way to yield the distribution of P. An 
outline of a diagrammatic proof of the cancellation of infra-red divergences 
is given in Appendix [A[ In Section ^ we examine the mean and covariance of 
P M at T = 0, since this is the base against which the T > results have to be 
compared. In Section || the corresponding T > quantities are computed. 
Various details of these calculations are relegated to Appendices [BfJB|. Our 
results and their implications are discussed in Section [5]. The reader who is 
not interested in the details of the computation may skip to this section. 



2 The Resummed Distribution 

In this paper we assume that a factorisation of the hard and soft parts of 
the cross section holds in the form given in eq. ( |1.2[ ), and then examine the 
soft part of the cross section, i.e., the distribution of P^. The generating 
function for all matrix elements needed for a resummation of the soft part of 
the one-loop order contribution is ||, 



Z(J,T) = J VA^ exp [-z J d A x Qi^F^ + J^A" 



(2.1) 



The integral over x$ in the exponent is carried out over a path in the complex 
time plane. The choice of this path is made so that the thermal averaging 
of the S-matrix elements is done properly. Detailed discussions of the path 
of integration can be found in ||. Some more information can be found in 
Appendix [A]. 

At this point it is useful to introduce some notation. For the scattering 
process, the initial and final momenta of the electron are p and p', and the 
initial photon has spacelike 4-momentum q (q 2 = —Q 2 ). For the annihila- 
tion process, the electron momentum is p, the positron momentum is p' and 
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the final photon momentum q' is timelike (q' 2 = Q 2 ). Details of the hard 
scattering kinematics are given in Appendix [B|. 

A proof of the factorisation implicit in eq. Q1.2Q can be obtained by con- 
sidering the leading perturbative corrections to the hard scattering process 
and taking the limit of the extra photon momentum, k, going to zero (see 
Appendix [A]) . For the processes we consider, this gives 

Here e is the positron charge. It is important to remember that this approx- 
imation is valid only in the infra-red. Hence, in order to use eqs. (|2.1|) and 
( |2.2|) consistently, we must truncate the mode expansion of the soft fields in 
the partition function. 

We perform this truncation in a Lorentz invariant fashion by writing the 
mode expansion — 

M$ = E [ d ~ k *W x (k)e lk - x + e,(k, X)a x (k)e~ lk - X ] , 

\L 1 (2-3) 

dk = 1 ^2n5(k 2 )e(k )e(riV 2 -V ■ k). 

Here e^fc, A) is the polarisation vector for a photon of momentum k and 
polarisation A, and a\(k) and a\(k) are the annihilation and creation opera- 
tors. V is a vector built out of the momenta of the hard particles which sets 
the scale of the hard process and C 1. The second theta- function in the 
Lorentz invariant phase space measure expresses the fact that the fields are 
soft. 

Note that such a factorised approach is useful only when there is a great 
difference between the hard and soft scales. In particular, we take all thermal 
fluctuations into the soft scale, and hence we must have Q 2 3> T 2 . For both 
the processes, we choose V to be the timelike vector q' (q' = p + p' in the 
hard-scattering kinematics). The arguments in this section do not depend 
on the choice of V. 

In the Feynman gauge, the generating function involves only Gaussian 
integrals. Performing them we obtain 

Z(J,T) = exp J ^J^k)G, v (k,u)r(-k)\ , (2.4) 
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where the thermal photon propagator is given by 



G^(k,u) = -g llu i [ -^--2m5(k 2 )B(k-u;T)Y (2.5) 

and the Bose distribution function is 

B(k ■ u; T) = e]k J T _ r (2.6) 

The Lorentz vector u (u 2 = 1) specifies the 4- velocity of the heat bath with 
respect to the observer [10|. In the rest frame of the heat bath, k ■ u = ko, 
and the Bose distribution takes its familiar form. This is a covariant way of 
expressing the fact that the presence of the heat bath selects out a preferred 
frame. 

The matrix element for n-photon emission is obtained by taking appro- 
priate derivatives of Z(J,T) with respect to the external current, giving 

M(h, •••,&„) = ZU^fle^J^h). (2.7) 

i=l 

Emission of photons corresponds to taking fc° > 0, and absorption to k® < 0. 
The squared matrix element, after summation over polarisations, gives a 
factor of J 2 for each photon. Finite temperature effects in the cross-section 
may be completely taken into account through an extension of the usual trick 
(see Appendix |A]) of writing a phase space element 

d<Pi = t^tt* (k 2 ) [9 (*°) + B{h ■ u- T)J Q( m ' 2 ~ Q ■ h). (2.8) 

The resummation now gives the probability density for radiation as 
dV 



d 4 P 



oo 1 „ n / n \ 

\zu n 2 y: - / n {dkwjuki)}* [p» - y,k )■ (2-9) 

n=0 n - i=0 \ i=0 / 



For further analysis we construct the Lorentz invariant characteristic function 
of this distribution — 

r dV 

S(a;) = jd*Pe iPx -^p = exp[R (x) + R T (x)]. (2.10) 
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The Fourier transformation decouples the delta function, allowing an expo- 
nentiation in the form shown. 

The logarithm of the characteristic function is defined by the expressions 

R = [ dkJ(k) ■ J(k){e ik x -l\ , 

J r ~ r , (2-11) 

R T = dkJ(k) ■ J(k)B(k-u;T)[e ik - x + e~ ik - x -2} . 



The subtractions come from \Z(J, T)| 2 and may be obtained from eq 
All infrared singularities cancel between the exponential terms and the con- 
stant piece in both the T = part, Ro and the T > contribution, R T . This 
is an extension of Weldon's proof of infrared finiteness in this model ||. 

The characteristic function, S, is a generating function for the moments 
of Pp. The logarithm, Ro + R T , is a generating function for the "connected 
parts" of these moments, i.e., the cumulants. The n-th cumulant of P^ is 
given by 



(2.12) 



x=0 



where the subscript c denotes a connected part. Each cumulant of P then 
becomes the sum of a T = and a T > part. 

Since the integral representation of R T is even in x, odd derivatives vanish. 
As a result, there is no T > contribution to the odd cumulants. The 
physics behind this is simple. At finite temperature the charged particle can 
remove momentum from the heat bath by absorbing a photon, as well as add 
momentum by emitting. In thermal equilibrium these two processes balance 
out, giving a distribution even in P. It is possible to pursue this further. 
The symmetry under x <-> — x comes from the phase space in eq. ( |2.S| ), in 
particular from the + B factor there. This is obtained by summing the 
weights for absorption, B(k )Q(— k ), and emission, [1 + B(k )]Q(k ). These 
two weights are related through the commutation relations of the photon 
field operators. In the ultimate analysis, this is the origin of the vanishing of 
odd cumulants in thermal equilibrium. 

It is useful to note the following result. Combining eqs. Q2.ll ) and (|2.12|) , 



we find that a cumulant of the resummed distribution is equal to the cor- 
responding moment of the distribution obtained from the soft part of the 
leading order in perturbation theory. 
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We write down explicitly the first cumulant, i.e., the mean radiated 
momentum — 

(Pa) = -t^-(R + R, 



(2.13) 

x=0 



(Pp) is given entirely by the T = piece, 

(P M ) = J dkJ(k) ■ J(k)kp (2.14) 

This is a Lorentz vector. We discuss it in detail later. 

We also write down the second cumulant (P^P^c = (P^Pv) — (P^)(Pv)- 
This is the covariance of the radiated 4-momentum — 



{P»Pu) e - - J" Q Xy (RO + Rt 



x=0 



V° +V T , where 



C = J dU(k) ■ J(k)k^ (2.15) 
V* = J dkJ(k) ■ J(k)2B(k ■ u; T)k^k u . 

The covariance is a symmetric Lorentz tensor of rank 2. The diagonal com- 
ponents of this tensor give the variance of the corresponding component of 
the radiated momentum. The off-diagonal elements specify the correlation 
between different momentum components. In particular, the space-space 
components quantify angular correlations in the radiated momentum. 

3 Zero Temperature 

The mean radiated momentum after resummation can be written in terms of 
all the vectors available in the problem, multiplied by appropriate coefficient 
functions (see Appendix 0). As a result — 



W = / dkJ x {k)j\k)k, = F x ^ + F 2 p . (3.1 




It is important to note that the frame dependence of all the tensors is taken 
entirely into the basis tensors. As a result, the coefficient functions are 
Lorentz invariant. One can obtain these and the tensor components quite 
trivially by contracting (P /t ) with each of the vectors on the right hand side. 
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For convenience, these Lorentz invariant contractions are evaluated in the 
Breit Frame (BF), i.e., the frame in which p and p' are equal and oppositely 
directed, 

h= =%vQm-m-n e ], 

h = i(P») = f vQ [(2* - 1)0! - fia] ■ 

Here a = e 2 /An and \I/ is defined in eq. (|B.1|) and the integrals Qi in eq. flC.7p . 
The final expressions are obtained in the BF and we have used the fact that 
Qi = Q't in this frame (see Table [I] in Appendix |D|). The equality of the two 
contractions translates into the equality of Fi and F 2 . In the hard-scattering 
limit (Q — > oo), we find 



Fi — F 2 — ^ 

7T 7 



21o,,(-)- 1 

m 



(3.3) 



where 7 is defined in eq. ( [B.4j ) . The logarithm of Q/m arises from that part 
of the phase space in which emitted or absorbed photons are collinear with 
the electron momenta. 

The components of the mean radiated momentum can now be easily eval- 
uated in any frame. In the BF, the only non-zero component is 

(P ) = 2 7 iq. (3.4) 

In an arbitrary frame, other non-zero components exist. The BF is a good 
frame in which to compare the T = and T > results because of the 
simplicity of the vacuum results in this frame. 

The variance of the radiated momentum has the tensor decomposition 

C = / dkJ x (k)j\k)k,k u = Y,F^Z; U , (3.5) 
J i=l 

where the four Lorentz tensors are listed in eq. ( |D.3| ) . The coefficient func- 
tions and tensor components can be computed as explained in Appendix 
0. 

We need to write the four contractions 



(3.6) 



J? = 


z 1 v° 

fJ,U ' [IV 


= 0, 








n = 


z 2 v° 

[IV ' [IV 


= f^ 2 Q 2 [2* fi 4 


- Oo- 


5 




Z 3 V° 

[IV ' [IV 


= § v 2 Q 2 


2^4 


-Oo- 




1° - 


z 4 v° 

[IV [IV 


= iyQ 2 


2O - 


§ 4 . 
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The final expression for each contraction is written in the BF. The coefficient 
functions are obtained by solving a linear system of equations, as explained 
in Appendix |D]. 

Since the matrix Z (Z,-,- = Z l „ 



terchange of Z 2 and Z 3 , 
scattering limit, we find 



Z l Z^, see eq. p.4| ) is symmetric under in- 
and 7° = one m ust have F® ~ ' 11 



Fn. In the hard 



71 

F° 



71 



7T 



vQ 2 



vQ 2 r 



log(§ 



- 1 



(3.7) 



Note that F2 = F3, as expected, and that collinear singularities appear in 
these two coefficient functions. 

The tensor components of the covariance can now be written down. We 
find that in the BF only the diagonal components survive. 



V° 
V° 
V° 



71 



V° 



71 



vQ 2 



2 



1 



-3 



(3.8) 



In a general frame, of course, there will be other components of the covari- 
ance. The coefficient functions, eq. ( |3.7| ), however, are Lorentz invariant, and 
may be extracted from measurements in any frame. 

Other important details to notice about these expressions are that there 
are no off-diagonal elements, i.e., there are no angular correlations. Also, 
the transverse asymmetry — V® 2 vanishes, and neither of the transverse 
components of the variance involves any large logarithms. Finally, the log- 
arithms in the other two components of the variance are equal, so that the 
longitudinal asymmetry Vq — V® 3 is free from these terms. 



4 Finite Temperature 

The finite temperature piece of the resummed distribution is even in P. As 
a result, the T > contribution to all odd cumulants vanish. In this section 
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we consider only the first non- vanishing cumulant, i.e., the covariance. It has 
the tensor decomposition 



V. 



dkJ x (k)J x (k)2B(k-u;T)k^k v = J2 F i% 



\IV1 



(4.1) 



where our choice of the seven Lorentz tensors is listed in eq. ( p.5|) . Recall 
that at T = only four tensors could be constructed. The three extra tensors 
are due to the extra vector u in the problem. The coefficient functions and 
tensor components can be computed as before. 

We begin by evaluating the contractions If = X^,V3,. These integrals 
are most conveniently performed in the rest frame of the heat bath (HBF). 
As discussed in Appendix |C| when Q ^> T the contractions are given to 
accuracy of order exp(— Q/T) by 



1 1 



1 7 



^-T 2 [fi 2 



^fT 2 [2 7 (fi x 



Q' 2 - 2#fi 3 ] , 
• 2^1) - (fi 



r 2 



ira 
T 1 



5 ' 

12 TT 

5 ' 



fi' 5 - 2VW 4 )] + 7 2 / ; 
2 7 '(fi' 1 + fi 6 - 2^i ) - (fi + fi 5 - 2^fi 4 )] + i 2 l 

2fi + (i ~ 2 7 ) Oi + (£ - 2 7 ') fi£ 



^T 2 [2^fi; 



fii — fig 



J 5 > 



^T 2 i [2*^ - fii - fie] - II 



(4.2) 



We have used the approximations given in eqs. ( p.8|) and ( |C.9| ). Note that 
the two pairs IJ and /J, and 1$ and If are related by interchange of primed 
and unprimed quantities. The angular integrals are listed in Table |I|. 

We write down only the leading behaviour (largest power of 7) of the 
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coefficient functions in the hard scattering limit- 



FT 
FT 



FT 



FT 
FT 



FT 



7raT 2 21og(l + A) 
— 3 7T ' 



naT z 1 



"1 — 2 72 



1 



iraT 2 1 

— q — 2 72 

<J 7 7+ A 



A(l -2A) + 2A 2 log(27) 

-(1 + 2A- A 2 )log7_ 
A(l - 2A) + 2A 2 log(27) 



;i + A) 2 log 7+ 



-(1 + A) 2 log 7 _ - (1 + 2 A - A 2 ) log 7+ 



iraT^jA-logjl + A)} 



A 



^[41og(27)-2], 
Zef- \ [2Alog(27) - A - (1 - A) log 7- - (1 + A) log 7+ ], 
^f- J[2Alog(27) - A - (1 + A) log 7 - - (1 - A) log 7+ ]. 



(4.3) 

In order to display these results in a compact form, we have introduced the 
notation 



x 



Ax cos x, 7± = 7„(l±a0, ^4 = 7+7- 



7 2 /3 2 sin 2 X , (4.4) 



where @ u , 7„ and x are defined in eq. flB.4|) . Since we are interested in 
the hard-scattering limit, we have only kept the leading power of 7 in these 
expressions. Note the appearance of log(27) = log(Q/m) in several of the 
coefficient functions. 

Although the integrals in eq. (|4.2|) were performed in the HBF, the co- 
efficient functions are invariant. However, the components of the covariance 
tensor are frame dependent. For ease of comparison with the T = results 
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m eq. 



V T 



V T 



V T 



v 33 



V T 



V T 



v 13 



we write the non-vanishing components in the BF- 

2 • 



f naT 2 



+ (1 + x) 2 log7_ 



+ (2A - 1 - 3x z 
-x) 2 log7+ , 




2-^log(l+A) 
A), 

1 + x 2 ) log 



|log(l 



[1 + A){1 



X 



(iraT 2 
f iraT 2 



f naT 2 



(1 



x 



+ (1 + x) 2 log7 
2j/? u sin x, 

2x log 



m z 
- + (1 



(2Ar 2 - 
2 log7_ 



-3) 



x 



[l + A)(l-x 



1 + x) 2 log7_ 



+ 2a;(A-2) + 
- (1 -x) 2 log7+ 



(1-x-) 



^-x/? u sinx- 



(4-5) 

A check is provided by the fact that the diagonal elements of this tensor are 
positive. 

We find several interesting features in eq. Q4.5I). Most obviously, there are 
some non-zero off-diagonal components of the covariance tensor. Recall that 
p, p' and u together define the reaction plane. Then all in-plane components 
of P are correlated with each other, but not with the off-plane component. 
Furthermore the transverse components of P no longer have equal variance. 
This is also due to the distinction between in-plane and off-plane components. 
Finally, all the logarithmically divergent terms arise only in the energy and 
longitudinal parts of the tensor. Exactly as at T = 0, the coefficients of the 
logarithms in the diagonal components are equal. 

An invariant comparison of the vacuum (T = 0) and T > results can be 
made using the coefficient functions. For this we need to convert the tensor 
expansion in eq. (|4.1|) into an expansion in terms of the set S 1 ^ (eq. |D.7| ) , since 
the first four tensors in this set are precisely equal to the tensor set used at 
T = 0. The coefficient functions in this basis are F? = J2j -PJS ,J . Using the 
explicit form for S given in eq. ( |JJ.8| ) , we see that Ff = F? for i — 1, • ■ • , 4. A 
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Figure 1: The regions of the photon phase space from which one obtains 
different logarithmic terms are included within dashed lines. The directions 
of the vectors in the problem are shown. Note that in the HBF the origin 
(marked with a cross) is not included in the integration because of infra-red 
cancellations. In the BF the ray parallel to u is removed. 



comparison of eq. ( |3.7|) with the first four expressions in eq. (|4.3j ) shows the 
change due to the temperature. Apart from the differences in the coefficients 
of the log(Q/m) terms, there are new log7± terms. Using the arguments in 
Appendix (see the discussion following eq. |C.9| ) it is easy to see that these 
are terms in \og(Q/T). 

The origin of the logarithms is shown in Figure [l]. The log(Q/m) terms 
arise from the same source at T = and T > 0. They come from the 
part of phase space where the photon momentum k is nearly parallel to p or 
p f . The log(Q/T) terms, on the other hand, arise entirely from the thermal 
distribution. This is clear from eq. (^4.1|) , since the j u dependence of the 
results can come only from the k ■ u term when the integrals are performed 
in the BF. Hence, these logarithms can be identified as due to the part of 
the phase space where k ■ u < T. Since this is a Lorentz invariant 0, we can 
transform this to the HBF, where it becomes equivalent to saying that the 
logarithms arise in the region of phase space where ko <T. 

The log(Q/T) terms never give large contributions. One way of under- 
standing this is to note that the actual singularity in the Bose distribution 
is removed through a cancellation between the real and virtual diagrams. 
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This is equivalent to removing the origin in the phase space of k. Thus, in 
the limit T — > 0, when the Bose distribution becomes concentrated on the 
origin, its contribution to the process actually vanishes. In other words, pos- 
itive powers of T always regulate the divergence in logT. Contrast this to 
the collinear parts. Along the rays parallel to p or p' in phase space there are 
actual singularities. These regions become available as m — > 0, and hence 
give rise to genuinely large contributions to various observables. 

When Q > T ^> m, then the covariances provide straightforward mea- 
surements of the temperature. Useful observables are the two asymmetries 

(P P ) C - (P 3 P 3 ) C {P\Pi) c - {P2P2U. ( . 

L ~ (P P ) C +(P 3 P 3 ) C ' T ~ (P 1 Pi) c + (P 2 P 2 ) c ' K ' } 



and the correlation functions 



C, u = - ^"' c =. (4.7) 
^(P,P,) C +(P U P„) C 

In the BF, A T , C i, C 03 and C\ 3 are precisely zero at T = 0, but become non- 
zero at finite temperatures. Each of these is proportional to (T/r]Q) 2 and 
hence need not be very small. As a result, they can yield measurements of the 
temperature. Furthermore, all other correlations remain zero even at finite 
temperature. The pattern of vanishing and non-vanishing correlations, or 
equivalently, off-diagonal components of the covariance tensor shows whether 
or not any observed effects are due to finite temperatures. As a result, these 
observations provide not only measurements of the temperature, but also 
cross checks that the effects do indeed come from thermal effects. 

We do not discuss numbers here since such estimates require matching 
of the perturbative and resummed parts of the cross section (TT|. It would 



nevertheless be interesting to carry out such a computation to see whether 
it is possible to design special purpose experiments to measure these effects. 



5 Discussion and Conclusions 

We have performed a resummation of soft photon emissions in two hard 
processes in finite temperature QED. These are the scattering of an electron 
by a hard spacelike photon, and the annihilation of an electron positron pair 
into a hard timelike photon. 
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The summation has been performed in a fully Lorentz-invariant man- 
ner. This was achieved by separating soft and hard photons invariantly, 
eq. ( |2.3| ), and by writing the thermal photon propagator invariantly, eq. ( |2.5| ). 
This second step reveals that finite temperature physics involves one extra 
Lorentz vector over the corresponding vacuum process |]10|| . This is a Lorentz- 
invariant way of stating that a heat-bath selects out a preferred frame. 

We showed that an invariant resummation of soft photon effects is possible 
and yields a perfectly finite distribution of the total radiated photon momen- 
tum, P. Thus Weldon's extension || of the Bloch-Nordsieck theorem || [12 



can itself be extended from energy distributions to the full 4-momentum dis- 
tribution. The net effect of finite temperatures is to produce correlations 
in the radiation pattern. By performing a tensor decomposition of the cu- 
mulants of P, we obtained a Lorentz-invariant comparison of the radiation 
pattern in vacuum and in a heat-bath in terms of certain scalar coefficient 
functions, eqs. ( |3.7| ) and (|4.3|). At finite temperature, large logarithms appear 
in unfamiliar places. 

The comparison of these results is simplest in the BF. In this frame there 
is a rotational symmetry around the common direction of the two fermion 
momenta at T = 0. Furthermore, there is a symmetry in the interchange of 
the two momenta. As a result, the mean momentum vanishes, leaving only 
a non-zero mean radiated energy, eq. ( |3.4| ). Also, there are no correlations 
between the different components of P, i.e., the covariance tensor is diago- 
nal. The transverse diagonal components of this tensor are equal due to the 
rotational symmetry. Large logarithms show up in the other two diagonal 
components, eq. (|3if ). 

At finite temperatures both these symmetries are spoilt due to the pres- 
ence of the extra vector u. The mean of P is not changed, since momentum 
can be absorbed from the heat bath as well as emitted. For a process in 
equilibrium, these two rates balance out, as seen from the symmetry of R T 
(eq. |2.1ip under x — > —x. However, the covariance tensor changes substan- 
tially (eq. |4.5| ). The extra vector, u, allows an invariant distinction between 
transverse momentum in the reaction plane and off it. The variances of these 
two components are no longer related by any symmetry, and hence can be 
different. All the in-plane components of P are correlated with each other, 
but not with the off-plane component. As a result, the covariance tensor is no 
longer diagonal. In particular, the covariance of the energy and longitudinal 
momenta contains a large logarithm, and may be easily measured. 
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We notice that the T > results contain two types of logarithms — either 
in Q/m or in Q/T. The source of the former is the same in the vacuum and 
finite temperature theories; they arise from parts of the phase space where the 
radiated (or absorbed) photon is nearly collinear with the fermion momenta. 
In the m — > limit photons can become precisely collinear with the fermions. 
This configuration gives a truly divergent contribution to various observables. 
As a result, the \og(Q/m) terms can become large. The log(Q/T) terms arise 
only at T > 0, and are due to soft parts of the phase space (in the HBF). 
The Born-level singularity at k = is removed by cancellation between real 
and virtual diagrams, and these logarithms never give large contributions. It 
appears that for hard scattering of thermal particles the log(Q/m) terms 
do not arise. In this sense, the scattering of external particles in a heat bath 
seems to be closer to the familiar situation at zero temperature. 

When including hard thermal loops, the formal considerations of section 
|2| change in two ways (see Appendix [A]). First, the phase space element in 
eq. (p.8|) has to be modified, by taking into account the absorptive part of 
the Braaten-Pisarski resummed photon propagator [13| . Second, the current 
in eq. (|2.2|) must be rederived using this resummed photon propagator in the 
leading perturbative correction. The formal results are otherwise unchanged 
provided the cancellation of infra-red singularities can be shown. 

Consequently, we expect the arguments in the rest of this paper to remain 
valid after Braaten-Pisarski resummation. In particular, odd cumulants will 
still vanish in the thermal part of the resummed distribution. There will 
still be seven coefficient functions for the covariance tensor, and the zero and 
non-zero components can still be identified by the methods developed here. 
Detailed computations including hard thermal loops are in progress, and will 
be reported later. 

The extension of these results to one-loop order QCD resummations is 
straightforward. We refer to the resummation procedure discussed in |l4j 



in 



the context of gluon radiative corrections to the Drell-Yan process. The phase 
space measure there is written in the Sudakov parametrisation. Nevertheless, 
it is easy to check that the resummation of the leading soft singularity in one- 
loop order corrections gives an exponentiation of J fJL {k)J fl {k). This current is 



very similiar to that in eq. ( |2.2|) . The only differences are that m 2 = in |14 
and that factors of a are replaced by a s Cp. Since no non-Abelian couplings 
appear in QCD to this order, all of Appendix [A] can be translated to that 
context with only these changes. 
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The upshot is that the resummation procedure outlined here also gener- 
alises to one-loop order resummation in finite temperature QCD. All our for- 
mula? can be carried over by the simple replacement a — > a s Cp. We expect 
that the logarithmic mass singularities can be absorbed into the structure 
functions of the external hadrons which undergo hard scattering inside the 
heat-bath. 

There are obvious phenomenological consequences of the computations 
presented here. Each off-diagonal component of the covariance tensor can 
be measured independently. A separate measurement of the temperature is 
obtained from each non-zero component. These estimates must agree with 
each other, and from that obtained by measuring A L and A T . A further 
experimental check is that all off-plane components of this tensor are indeed 
zero. These conditions being fulfilled would verify that the observations are 
due to thermal effects, and simultaneously give a measurement of the tem- 
perature. In summary, soft photon resummations yield a new thermometer 
in high energy physics. 
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A The Resummation Procedure 

Here we give a diagrammatic proof of the resummed result to order a for 
the processes considered in the text using the techniques of 0. The hard 
scattering cross sections with soft thermal photon contributions can be ob- 
tained in perturbation theory. The hard particles are external and do not 
thermalise in the heat bath. The interaction Lagrangian is 



where u} is written in terms of the hard part of the photon field and C} in 
terms of the soft part — 



C I = C h I + £ 



s 



(A.l) 



C) = j» • Aftx) and L\ = j» • A*{x). 



(A.2) 
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In this appendix we use the notation ■ = J d xj ,J, (x)A tJt (x) and the 
current j M = ieip^^ip. In order to make the computation transparent, we 
quantise the fields in a box of volume V and take infinite volume limit at the 
end. Thus, 

A K X ) = 12 v /^ry [ a i e M exp(-i/c; ■ x) + aje*(/c z ) exp(i^ ■ x)\. (A.3) 

The prime over the summation denotes that it involves only soft modes. The 
continuum Lorentz invariant version of this is given in eq. ( |2.3| ). 

In terms of the transition matrix, T, and the density matrix for the 
thermal fields, p, the inclusive cross section is 

a = ^Tr(pTTt) 

= Z^F exp(-^/3n^) ({m},e'|T|{n},e,7*) c 

{m},{n},e' I 

where Z = Hi[{ — ex p(~ and iT = exp[i J dxCi(x)} — 1. 

(A.4) 

The subscript c in the matrix element denotes the connected part. The label 
{n} stands for the set of occupation numbers for each mode of the soft fields, 
and F is a flux factor. Note, p = exp(—/3H)/Z where Z is the full partition 
function of the interacting system. The process of taking connected parts in 
the second line absorbs the interaction effects and leaves only a normalisation 
by the ideal gas partition function Z . The density matrix for electrons in 
the heat bath is disregarded since the hard particles are not thermal averaged 
and because there are no electron loops to lowest order. 

The soft photon contribution to the hard scattering cross section can be 
obtained from the following part of the transition matrix — 



(A.5) 



where T denotes time ordering. This expansion is made keeping the order a 
term in the hard scattering and all orders in the soft thermal photons. After 
Wick contractions, this gives two types of graphs — real soft photon emission 
and absorption, and virtual soft photon corrections. For every field which is 
thermal averaged the contraction gives a thermal propagator. 
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p 



Figure 2: The virtual corrections discussed in eq. ( A.7Q . The first diagram 



corresponds to the lowest order process in a v , and the other three to the 
parts in M v . 

To order a the virtual soft photon contribution (Figure to the hard 
scattering cross section can be obtained from 

T v = l -T{f-Alf.A s a jP.A s p ), (A.6) 

by Wick contracting only the soft thermal photons. To this order there 
are no thermal fermion loops. Since there are no external soft photons in 
the virtual part, only the diagonal elements of the transition matrix will 
contribute. Taking the continuum limit, we obtain the result in an arbitrary 
frame as 

a v = * J ^0^^{P 2 - m>)&( P > )(2nm P + q - p') 

xu( P y (q)u(p)]*[u(p')M v u(p)], 
where M v = e^(k,p,p') + Z(k,p')j±^j(q) + ^ q )j±-E{k,p). 

(A.7) 

Here S(/c,p) and T^(k,p,p f ) are the usual self-energy and vertex functions 
calculated at T = 0, only with the photon propagator replaced by the thermal 
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Figure 3: The real photon emission and absorption corrections of a r in 
eq. ( A.lUp . Absorbed soft photons come in from the left, emitted photons 
exit to the right. Note that the distinction between emission and absorption 
is included into the phase space measure. 



photon propagator of eq. ( |2.5p . The electron propagators do not have thermal 
parts, since they are assumed external. In short, this expression is identical 
to the T = result but for the thermal photon propagator. This remains 
true even if one considers the Braaten-Pisarski resummed hot thermal loops 
[ PHI , although the form of the propagator then changes. 

The real soft photon emission and absorption graphs (Figure ||) upto the 
same order are obtained from the term 

T r = T(j^A h ^r-A s a ). (A.8) 

The matrix element of this operator is proportional to ({m}|v4* |{n}}. Only 
the off-diagonal elements of the transition matrix contribute to this cross 
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section. Using the identity 



1 



Zi 



(A.9) 



hi 



(where (3 = 1/T) we get the following cross section for real photon emission 
and absorption in the 3 limit — 



e_ 
F 



d 4 p' 
(2vr) 4 



27r5{p' 2 - m 2 )Q{p' ) d<f>(k) 



x (2ir) 4 5 4 (p + q—p' — k) 



where M r = j(q)-rr 



1 



u{p')M r u{p) 



(A.10) 



p 1 + 1/t — m 



m + m 



rn 



and the phase space element for real thermal photons is given by eq. (|2.8| )f[ 
This phase space for T > changes because of the additional absorption 
processes which do not exist in the vacuum theory Again, the electron 
propagators do not have thermal parts. The important point is that in both 
the real and virtual corrections, the matrix elements are the same as at 
zero temperature. Finite temperature effects are subsumed into the thermal 
photon phase space and propagator. 

We note that the squared matrix element appearing in a r amounts to 
taking the one loop corrections to the basic process with a cut on the soft 
photon line. The phase space thus comes from the absorptive part of the 
thermal photon propagator. Upon performing Braaten-Pisarski resumma- 
tion, the real photon phase space element will then be changed in obvious 
ways. 

The soft photon limit is obtained by setting k to zero in the numerator 
of the matrix elements in a v and a r . In this limit we get 



[u(p')fiu 



[uip'WvU 



h(k)nkw 



u(p')M r u{p) 



u(p')£(q)u(p) 



(A.ll) 



where J M is the current in eq. ( |2.2|) . This is the origin of the factorisation in 
eq. (|1.2|) . Taking higher order terms from eq. ( |A.5| ) gives the exponentiation 



section ^ a Lorentz invariant cutoff on the soft modes has also been included with 
this. The rest of the phase space element was also written down in H. 
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of the order a result. It is possible to show that a factorisation in the form 
shown in eq. ([A. 11 ) holds with higher powers of (e^J^) 2 . As a result, we 



obtain the squares of the matrix elements shown in eq. ( j2.7p , and hence the 
generating function of eq. (|2.1|). A full proof will be given elsewhere. 



B Hard-Scattering Kinematics 

The problem considered here contains three Lorentz vectors, the two fermion 
momenta, p and p', and the velocity vector of the heat-bath, u. With three 
vectors we can construct three Lorentz invariants. Since p 2 = p' 2 = m 2 and 
u 2 = 1, we can work with the dimensionless invariants 

1 I ; l 

— p-u = 7, — p 1 ■ u = 7', and — -p ■ p' = \F (B.l) 

m m m 2 

For T = 0, the problem does not involve u and the only invariant available 
is ^. 

For the electron scattering problem, it is a convention to use the square 
of the vector q = p' — p as one of the invariants. This is a measure of the 
scale of the process. We have the relation 

Q 2 = - q 2 = 2m 2 (^-l). (B.2) 

The limit of a hard scattering, Q 2 — > 00, corresponds to taking \1> — > 00. 
The momentum q' = p + p' is closely related. In the annihilation process this 
corresponds to the final photon momentum, and its square to its invariant 
mass. We have, 

q' 2 = 2m 2 (V + l)^Q 2 . (B.3) 

Note that q is space-like and q' time-like, and q.q' = 0. We choose to use q' 
to set the scale of the hard scattering process. 

We shall find it most convenient to work in one of two frames, either the 
rest frame of the heat-bath (HBF) or the Breit frame (BF). In this latter 
frame, p and p' have equal and opposite 3-momenta, and their common 
direction conventionally defines the z-direction. For the annihilation process 
the BF is just the rest frame of the final off-shell photon. The reaction 
plane in the HBF is chosen to be the xz-plane. Then this is also the plane 
containing p, p' and u in the BF. 
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We find it useful to parametrise different momenta in the BF as displayed 
below 

p = m7(l,0,0j), 

P > = 7717(1,0,0,-^), 

u = 7 u (l,/5 u sinx,0,/3 u cosx), 
g' = 2m7(l, 0,0,0). 

Note that the BF is the rest frame of q' . The parametrisations in the HBF 
can be obtained by boosting these so as to obtain u = (1, 0, 0, 0). 

The quantities 7, 7^ and \ are Lorentz invariant. It is easy to see that 



- 1 /* + 1 

" = vW and 1 = v— • (a5) 

As a result, the hard-scattering limit corresponds to taking 7 - ► 00. In 
addition, we can write an explicitly Lorentz-invariant form for j u , 

7„ = "^£±£1. (B.6) 
An explicitly invariant form of x can a ls° be written down. 



C The Integrals 

In this appendix we investigate integrals of the form 

I = J j^y 4 2Tr6(k 2 )&(k )Q ( m ' 2 -q'-k) 2B(k ■ «; T)F(k), (C.l) 

where B is the Bose distribution function and F(k) is a scalar involving 
(possibly) p, p' and u. The second O-function on the right restricts the 
phase space to soft modes. Since we are interested in the region 77 <C 1, we 
shall freely use the hard scattering kinematics of Appendix 0. 

In the BF this phase space constraint takes the particularly simple form 

© (V 2 - 4 ■ k) = ®{2m^-k Q ). (C.2) 
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The only constraint is on the values of ko- For T = 0, the Bose distribution 
drops out, and the only possible angular dependences are in F{k). Hence, 
all the relevant integrals are best done in the BF. 

For T > 0, the Bose distribution picks up angular factors in the BF. 
Hence, it is simpler to perform these integrals in the HBF. However, the 
constraint now cuts off part of the angular region in phase space. Fortunately, 
the analysis of this problem is not complicated. The constraint can be written 
as 

e (v 2 - 4 -k) = e (^1 - mi - /W • k)) , (c.3) 

where boldfaced symbols with carets denote unit 3-vectors. 

The quantity within brackets on the right hand side lies between 1 — \(3 U \ 
and 1 + \p u \. As a result, the constraint picks a non-trivial angular region of 
integration only for k_ < k < k + , where 



= Imrf 



k + = 2mr]'- 



l + |Ax| 



(C.4) 



For k < the full angular region is allowed, and for k > k + , no phase 
space volume is allowed. 

Now, the integral / can be split up into an integration over < ko < k- 
(Ii) and another over k_ < k < k + (I 2 ). The angular integral in the first 
part is over the full sphere k, whereas that in the second part is restricted 
by the constraint. 

This second integral becomes 

1 r k + 



8vr 3 



/ + kdkB(k;T)f{k), (C.5) 

Jk- 



where f(k) is obtained from F(k) by the appropriate angular integration. 
Since B(k;T) is monotonically decreasing, we have the bound 

\h\ < B(k 8 ~; T) [ + kdk\f{k)\. (C.6) 

For all the integrands we are interested in, the growth with 7 of the integral 
on the right hand side is polynomially bounded. Since B{k^;T) decreases 
exponentially with 7, I2 is exponentially small in the hard scattering limit. 
As a result, it can be neglected in comparison with Ii. 
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The angular integrals in the first part are unconstrained. Hence we find 
it useful to define the integrals 

a 



J<KlFi(k), where 60.= *™*^ , (C.7) 

with the convention that Fi(k) is rendered dimensionless by dividing each 
factor of p or p' by m and k by fco- In the full computation we need many 
forms of Fi(k). These, and the resulting angular integrals are listed in Table 

B 

For Q — > oo, we can now write 

I w ( k ~ kdkB(k;T)F(k,9,(P). (C.8) 

2ir 2 Jo 

In actual applications, we shall choose k- to be large enough that 

dkkB(k-T) w . (C.9) 

6 

A numerical evaluation of the incomplete Bose integral on the left shows 
that this approximation is valid to an accuracy of one part in thousand for 
k_ pa 8.6T. Now if we were to have T <C Q, then we must either take 
i] <C 1 or 1 — \/3 u \ <C 1 in order to satisfy this condition. Since rj is some 
fixed number, for Q/T ^> 1, it is possible to satisfy the constraint even for 

lu > 1. 

In short, in the hard-scattering limit, we can approximate all integrals of 
the form in eq. ( |C.1| ) with exponential accuracy to obtain 

h = ^p. (CIO) 

The angular integrals Qi are frame dependent, and are listed in Table [I]. 



D Tensor Decompositions 

The cumulants are symmetric Lorentz tensors which can be written as an 
integral I^^-^ over a vector k, containing some vectors Vi, V 2 , etc, in the 
integrand. As a result, they admit expansions of the form — 
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Integral 


Fi 


HBF 


BF 


^0 


1 


1 


1 




[mk /p ■ k] 




£(0)/7 


n 2 


[m 2 kl/{p-kf\ 


i 




n 3 


\m 2 k} i lp ■ kp' ■ k] 


HP*)/* 




n A 


\p' • k/p ■ k] 


[0 + $L(/3)/ 7 ]/(/3 7 ) 2 


2L(0) - 1 






[e + AL(/?)]/(/3 7 ) 4 , 




^5 


\p' ■ k/p ■ k} 2 


f = $ 2 + 2(f) 2 - r 2 , 
A = r 2 -0(0-2$/ 7 ). 


4 7 2 + 1 - 4L(/3) 


^6 


mko\p' ■ k/(p ■ k) 2 } 


[$ + 0LO3)/ 7 ]/(/3 7 ) 2 


2 7 -L(/3)/ 7 



Table 1: Definitions and values of the angular integrals in the HBF and BF. 
The integration is over the full sphere k.We have used the mnemonics (a) 
L(x) = log((l+x)/(l-x))/2x, (b) p 9 = yfl- l/^ 2 , (c) = 7 y-*, (d) $ = 
\P 7 — 7 ' (e) <£>' = ^ 7 ' — 7 and (f) T = /3/3' 77 '. are defined by interchanging 
p and p'. In the HBF, their values are obtained by interchanging primed and 
unprimed quantities. In the BF, f2j = Q[, since each integral is symmetric 
under (3 — > —/3. 
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The set of basis tensors Q l consists of all possible symmetric rank-n tensors 
which may be written using the vectors Vi, V2, etc, and the metric tensor 
g^ u . We choose the basis tensors to be dimensionless. 

The coefficients Fi are Lorentz scalars, and given by the solution of the 
set of linear equations 



= ^QijFji where X { 



Ql TfJ-l-'-fln 



Q 



1.1 



(D.2) 

The matrix of contractions of the basis tensors, Q, is symmetric. The inte- 
grals Xi are Lorentz invariant. 

The only vectors involved in the problem at T = are p and p'. The 
tensors necessary for the mean momentum are then p/m and p'/m. For the 
variance, we need to consider the four tensors 



Z 1 



Z 



[IV 



Z 4 



■72 _ ViiVv 

P^Pu + VvV'a 



(D.3) 



The corresponding matrix of contractions is 

f 4 1 1 
1 



1 
1 



^2 
1 
1 



1 
1 
1 

1 + y 2 

J 



(D.4) 



The linear system to be solved is ZF = I. 

For the covariance at T > 0, the basis tensors must be built out of three 
vectors, p, p' and u. It is more convenient to work in terms of u and the 
two vectors p = p — (p ■ u)u and p' — p' — {p' ■ u)u. These two vectors are 
orthogonal to u |T(J. There are seven basis tensors, and they can be chosen 
to be 



r 1 

fJ,U 



T 3 
T 5 

fJ,U 

T 7 

[1U 



P~'j>'v 
2u ■ p 



T 2 



T 4 



t: 



PuPu 

2p-p' 
2u ■ p 



(D.5) 
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With this choice of tensors, the matrix of contractions is 



3 


-7 2 /3 2 




_ SL 











-7 2 /? 2 


7 4 /? 4 


2 











- 7 ' 2 /?' 2 


2 


















7 2 /5 2 


7'W 


r 2 + 2 















\|/ 


2^ 2 














1 























p 2 






















2 7 y 





(D.6) 



V 

We have used the notation for and V as given in the caption of Table [I]. 

For comparison with T = 0, it is more convenient to express the results 
in terms of the tensor set 



S 



fJbV 



Z l „„ for i 



S e . 



l,..-4, 



2u ■ p : 



2u ■ p 



(D.7) 



These can be expanded in terms of the set T\ We have T^ u = S^S*^, where 



' hi/ i 



/ 


1 











-1 








\ 







1 








7 2 


-2 7 2 















1 





7' 2 





-2 7 /2 















1 




77' 


77' 






























1 
























-1 


1 







\ 














-1 





1 


/ 



(D.8) 
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